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Abstract 

In this paper we investigate the local and global spectral properties 
of the triaxial rigid rotator. We demonstrate that, for a fixed value 
of the total angular momentum, the energy spectrum can be divided 
into two sets of energy levels, whose classical analog are librational 
and rotational motions. By using diagonalization, semiclassical and 
algebric methods, we show that the energy levels follow the anomalous 
spectral statistics of the one-dimensional harmonic oscillator. 

1 Introduction 

Berry and Tabor 1,2 have shown that for a generic quantum system, which is 
classically integrable, the spectral statistics of energy levels follow the Pois- 
son ensemble, while for a system whose classical analog is chaotic the spectral 
statistics follow the Gaussian Orthogonal Ensemble (GOE) of Random Ma- 
trix Theory (RMT) (for a review see Ref. 3). Nevertheless, there are some 
exceptions to this classical-quantum correspondence: the best known case is 
the harmonic oscillator one. 4 ' 5 Recently we have discussed another patholog- 
ical case: the classically integrable triaxial rigid rotator. 6 -* We have numer- 
ically computed the energy levels and then calculated the nearest-neighbor 
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spacing distribution P(s) and the spectral rigidity A^(L). We have found 
that -P(s) shows a sharp peak at s = 1 while A 3 (L) for small values of L 
follows the Poissonian predictions and asymptotically it shows large fluctu- 
ations around its mean value. 6 Such behavior of spectral statistics is quite 
similar to that of a one-dimensional harmonic oscillator. In this paper we 
have calculated the exact energy levels, the semiclassical ones and those ob- 
tained by an algebraic harmonic approximation. We show that, for a fixed 
value of the total angular momentum, the Hamiltonian of the triaxial top 
describes a sort of nonlinear pendulum and its semiclassical energy spectrum 
can be divided into two sets of energy levels, corresponding to librational 
and rotational motion respectively. These energy levels are close to the ex- 
act ones apart near the separatrix. Around the minima and the maxima 
of the classical energy, the exact energy levels are described accurately by 
approximate algebraic formulas. 

2 Triaxial Rigid Rotator 

The classical Hamiltonian H of the triaxial rigid rotator is given by 

H = \ (aJ{ + bJ 2 2 + cJt) , (1) 

where J = ( J±, J 2 , J3) is the angular momentum of the rotation and a = l/Ii , 
b — 1/ 12, c — 1/ 1 3 are three parameters such that Ii , /g and I3 are the princi- 
pal momenta of inertia of the top. 7 ' 8 Here we choose a < b < c. The quantum 
Hamiltonian H is obtained by replacing components J m of the angular mo- 
mentum, in the classical expression (1), by the corresponding quantum op- 
erators Ji, J 2 and J 3 . These obey the commutators [Jf, = i£e,j,mJm an d 
commute with the Casimir operator J 2 = J\ + J| + J|, namely [J 2 , Jk] = 0. 
In view of the latter equation, a standard choice to deal with the angular- 
momentum dynamics consists in diagonalizing the two operators J 2 and 
J 3 within the basis {\j, k) : \k\ < j} such that Jz\j,k) = kh\j,k) and 
J 2 \ji k) = j(j + l)U 2 \j, k), where the index j labels the angular-momentum 
representation. After recalling that J±\j,k) = [(j =F k)(j ± k + \)] l / 2 h\j, k± 1) 
(with J± = Ji ± i J 2 ), the non-zero matrix elements of the quantum Hamil- 
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tonian H in the basis \ j, k) are given by 

(j, k\H\j, k) = ^(a + b)(j(j + 1) - k 2 ) + ^ck 2 , (2) 

(j,k\H\j,k + 2) = {j,k + 2\H\j,k) = 

= y (a - b)y/(j -k)(j-k- l)(j + k + l)(j + k + 2) . (3) 

For a fixed value of j, the Hamiltonian matrix can be decomposed in four 
submatrices, corresponding to different classes of symmetry. 9 The numerical 
diagonalization of the matrix gives the "exact" energy levels of triaxial rigid 
rotator. 6 



3 Semiclassical description 



It is important to observe that, also if the "exact" energy spectrum of the 
triaxial rigid rotator cannot be analytically determined, the classical Hamil- 
tonian of the triaxial top is integrable. 7 ' 8 We use the classical integrability of 
our triaxial top to calculate the semiclassical energy spectrum via torus quan- 
tization. The first step is to find out the action variables of the Hamiltonian. 
By using the Deprit transformation 8 (this is also known, at the quantum 
level, as Villain's representation of angular momentum 10 



Ji = y/j 2 - Jf sin 9 , J 2 = yjj 2 - Jf cos , (4) 
where J 2 = J\ + J 2 + J 2 , the classical Hamiltonian can be written as 

H = \aJ 2 + \ Uc-a) - (6 -a) cos 2 ^1 J 3 2 + -(& - a) J 2 cos 2 6 , (5) 

where the dynamical variable J3 is the action variable of the angle 8. 8 By 
performing another canonical transformation the angular dependence of the 
Hamiltonian H = H(J, J3, 0) can be removed. Thus the Hamiltonian can be 
written as a function of two action variables: the total angular momentum 
J and a new action variable I given by 



j _ J [° 2 jr-a-(b-a) cos 2 6 
71 J61 \ (c — a) — (b — a) cos 2 6 



d6 



(6) 
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The Hamiltonian H is now an implicit function of two action variables: H = 
H(J, I). It means that the system is integrable and the classical trajectories 
of the 6-dimensional phase-space are restricted to a 2-dimensional torus. 7 For 
a fixed value of J the Hamiltonian (5) describes a sort of nonlinear pendulum 
where the extrema of integration are 



Q\ = arcos 
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for < H < (6/2)J 2 and it corresponds to librational motion, while they are 

X = , 6 2 = n , (8) 

for (6/2) J 2 < H < (c/2)J 2 and it corresponds to rotational motion. Note 
that H = (6/2) J 2 is the energy of the separatrix between librational and 
rotational motions. The semiclassical (torus) quantization 11 of the energy is 
performed by setting 

j + l)h, i = kh, (9) 



J 



where j and k are integer quantum numbers such that k = — j, — j + 1, j — 
1, j and j > 0. Note that only recently the corrections to this torus quan- 
tization, namely the higher order terms of the WKB series of the angular 
momentum, have been investigated. 12 For a fixed value of j there is a set of 
2j + 1 energy levels but the semiclassical energy spectrum E^ k , given by the 
implicit equation 
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(10) 



where ef h := E s A/[(j + l/2) 2 /i 2 ], has j couples of degenerate levels because 



Ef_ k . The lowest energy level of the set is = ati 2 (j + l/2) 2 /2 
while the higher is Ef^ = ch 2 (j + l/2) 2 /2. 



rpsc 



4 Numerical results: spectral fluctuations 

In Figure 1 we compare the "exact" energy spectrum, obtained with the nu- 
merical diagonalization 6 of the matrix Hamiltonian given by Eq. (2) and Eq. 
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(3), with the semiclassical one, calculated using Eq. (10). As expected, the 
semiclassical quantization cannot take into account the broken degeneracy 
of pairs of "exact" energy levels. Moreover we find that the semiclassical en- 
ergy levels are more accurate in the higher part of the spectrum (rotational 
levels) . 

Berry and Tabor 1 ' 2 have shown that for classically integrable systems the 
spectral statistics P(s) and A 3 (L) are expected to follow Poisson predictions: 
P(s) = exp (— s) and A 3 (L) = y|. These predictions are based on the analysis 
of the semiclassical energy spectrum of a generic integrable Hamiltonian with 
more than one action variable. Thus, to test the Berry- Tabor theory, is much 
more appropriate to investigate spectral statistics of the semiclassical energy 
levels than spectral statistics of "exact" energy levels, as done in Ref. 6. 

In Figure 2 we plot P(s) and A 3 (L) obtained from the semiclassical energy 
spectrum. The level spectrum is mapped into unfolded levels with quasi- 
uniform level density by using a standard procedure described in Ref. 13. 
P(s) is the distribution of nearest-neighbor spacings Sj = (-E«+i — Ei) of the 
unfolded levels Ei. It is obtained by accumulating the number of spacings 
that lie within the bin (s, s + As) and then normalizing P(s) to unit. The 
statistic A 3 (L) is defined, for a fixed interval (— L/2, L/2), as the least-square 
deviation of the staircase function N(E) from the best straight line fitting it, 
where N(E) is the number of levels between E and zero for positive energy, 
between — E and zero for negative energy. Figure 3 shows that the spectral 
distribution P(s) has a peak near s=l and nothing elsewhere. The spectral 
rigidity A 3 (L) follows the Poisson prediction A 3 (L) = L/15 for small L but 
for larger values of L it gets a constant mean value with fluctuations around 
this mean value. 

The results of Figure 2 and Figure 3 can be explained in the following way. 
As shown in the previous section, fixing the angular momentum J the triaxial 
rigid rotator is one-dimensional: in the case of systems with only one action 
variable the energy spectrum is locally harmonic and, after unfolding, the 
spectral fluctuations have the pathological behavior of the one-dimensional 
harmonic oscillator, thus the spacings are all close to one. Moreover, it is 
important to observe that both "exact" and semiclassical levels have a linear 
grow apart near the separatrix where they cluster. This effect is shown in 
Figure 3 and Figure 4 where we plot the true nearest neighbor spacings 
(without unfolding) of semiclassical and "exact" energy levels as a function 
of the energy for a fixed value of the angular momentum. The spacings of the 
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semiclassical energy levels decrease quasi-linearly with the energy up to the 
separatix and then they increase quasi-linearly The same behavior is found 
for the spacings of the "exact" energy levels. In particular, in Figure 4 we 
plot the spacings of the exact energy levels corresponding to the four classes 
of symmetry of the system: the total Hamiltonian matrix is decomposed in 
the direct product of 4 submatrices by considering the parity of the quantum 
number k: even (E) or odd (0), and the symmetry of the state: symmetric 
(S) or anti-symmetric (A). So the submatrices are labelled as follow: (E,S), 
(E,A), (0,S), (0,A). These are the classes of symmetry of the system. 6 ' 9 Note 
that the spacings of each class of symmetry have the same trend but there are 
some fluctuations, in particular around the extrema of the energy spectrum. 
We have verified that these results do not depend on the choice of the inertia 
parameters of the asymmetric rigid rotator. 

5 Algebraic Harmonic Approximation close 
to minima and maxima: single levels 

The investigation of the triaxial-top classical Hamiltonian with a < b < c re- 
veals that the minimum-energy configuration is doubly realized by ( J±, J 2 , J3) ~- 
(±J, 0,0), where J is the radius of the sphere Jf + Jf + J| = J 2 and 
represents a constant of motion. On the other hand, the configurations 
(Ji, J 2 , J3) = (0,0, ±J) and (J 1; J 2 , J 3 ) = (0, ±J, 0) are recognized to be 
maxima and saddle points, respectively. By using the Casimir operator J 2 , 
we notice that the quantum Hamiltonian H can be rewritten in the two forms 



a.J 2 + (b- a)J 2 + (c- a)Ji\ , (11) 
cJ 2 -(c-a)Jl-(c-b)Jl\ , (12) 



close to the minima and the maxima, respectively. The observation that, 
classically, the energy minimum is reached for J 2 = J3 = 0, and J\ = ±J, 
leads to approximate Ji as Ji ~ ±J[1 — (Jf + Jf)/2J 2 ]. This implies, at 
the quantum level, the inequalities (J s ) <^ J [where the radius J is given 
quantally by J 2 = j(j+l)h] concerning the expectation values of J s (s = 2, 3) 
and makes natural to adopt the perturbative scheme (see Ref. 14) in which 
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the operator J 2 and J 3 can be treated as canonically conjugate variables. 
Such a feature follows from 



[J 2 , J 3 ] ~ iJ , 



(13) 



wheras 



[i 3 , A] = iJ{j 2 /J) ^ , [4 i 2 ] = * j( j 3 / j) - , 



(14) 



are considered vanishing when compared with (13) since ( J s ) / J << 1. Then, 
after labeling by the integer k the eigenstates of the harmonic oscillator term 



(c-a)[Ji 



■W 2 J 2 ] contained in H, where W 2 



(b — a)/(c— a), its eigenvalues 



are found to be JW(2k + 1) so that the spectrum of H reads 



Tphar 



aj(j + l)h z + Jj(j + l)hJ(b - a)(c - a) (2k + l)ft 



(15) 



The approximation just implemented is not consistent with the spectrum 
expected for c — > a. In fact, in such a case the Hamiltonian becomes H = 
[cJ 2 + (b — a) J 2 }/ 2 thus entailing however a quadratic spectrum E = [cJ 2 + 
[b — a)k 2 h ]/2. Similarly, the approximation for the maxima is obtained by 
noticing that J 3 ~ ±J[1 - (J 2 + J|)/2J 2 ], and 



[Ji, J 2 ] ~ iJ, 
[J 3 , A] = tJ(J 2 /J) ~ , [J 2 , J 3 ] = iJ(Ji/J) ~ . 



(16) 
(17) 



In this case, the harmonic oscillator term (c — &)[J| + fi 2 Ji], with f2 2 = 
(c — a)/(c — b), contained in H leads to the spectrum 



rphar 



cj(j + l)h 2 - y/j(J + l)hy/(c-a)(c-b) {2k + 



(18) 



which differs from Eq. (15) owing to the interlevel separation. In Figure 5 
we compare "exact" energy levels with those obtained using the harmonic 
approximation. The results show that indeed the harmonic approximation is 
quite accurate, also better than the semiclassical one, close to the minimum 
and the maximum of the energy. 
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Conclusions 



In this paper we have investigated the energy spectrum of the triaxial rigid 
rotator, a very important model for both atomic and nuclear physics. We 
have found that its semiclassical energy spectrum is spanned by two quantum 
numbers: the quantum number j of the total angular momentum J and the 
quantum number k of the other action variable / of the system. For a fixed 
value of j the semiclassical energy spectrum depends on the quantum number 
k and can be divided into librational and rotational energy levels. Moreover, 
the spacings of both "exact" and semiclassical energy spectra show a quasi- 
linear behavior as a function of energy: nearest neighbour spacings have a 
quasi-linear decrease up to the energy of the separatix and then a quasi- 
linear increase. We have also shown that the semiclassical energy spectrum 
is in good agreement with the "exact" one, but cannot take into account the 
broken degeneracies of pairs of "exact" energy levels. From our analysis it 
follows that the spectral fluctuations of the asymmetric top have the patho- 
logical behavior of the one-dimensional harmonic oscillator. This effect has 
been verified by calculating the spectral statistics P(s) and A 3 (L). Finally, 
we have deduced analytical formulas, based on an algebric approximation of 
the quantum Hamiltonian, which describe remarkably well the energy levels 
of the triaxial rotor close to extrema of the energy. 

V.R.M. is greately indebted to M.V. Berry, S. Graffi and A. Turchetti for 
many suggestions at different stages of the paper. L.S. thanks T. Prosen for 
enlightening discussions. 
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Fig. 1: "Exact" (left) vs semiclassical (right) energy spectrum of the triaxial 
rigid rotator with a fixed value of the angular momentum j. Parameters: 
a = 1, b = v^2, c = a/3 and k — 1. Energy in arbitrary units. 
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Fig. 2: Nearest neighbor spacing distribution P(s) and spectral rigidity 
As(L) of the semiclassical energy spectrum of the triaxial rigid rotator with 
a fixed value of the angular momentum j. The dashed lines are the Poisson 
predictions: P(s) = exp (— s) and A 3 (L) = L/15. Parameters as in Figure 1 
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Fig. 3: True spacings (without unfolding) as a function of energy. Semiclas- 
sical energy levels of the triaxial rigid rotator for two values of the angular 
momentum j. Parameters and units as in Figure 1. 
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Fig. 4: True spacings (without unfolding) as a function of energy. "Exact" 
energy levels of the triaxial rigid rotator for a fixed value of the angular 
momentum j. Four classes of symmetry: square for class (E,S), circle for 
class (E,A), diamond for class (0,S), triangle for class (0,S). Parameters 
and units as in Figure 1. 
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Fig. 5: "Exact" (left), semiclassical (middle) and harmonic-approximation 
(right) energy spectrum of the triaxial rigid rotator with a fixed value of the 
angular momentum j. Parameters and units as in Figure 1. 
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